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ABSTRACT 

We have discovered 10 periodic X-ray sources from the 1 Ms Chandra ACIS observation of the Limiting 
Window (LW), a low extinction region (Ay ~ 3.9) at 1.4° south of the Galactic center. The LW provides a 
rare opportunity of studying the Galactic Bulge sources without obscuration from molecular clouds. Using the 
early 100 ks Chandra exposure and the HST observation that covers five of these sources, we have reported 
three of them as candidate CVs o r accreting binaries, based on their blue optical colors, excess Ha fluxes, and 
high X-ray-to-optical flux ratios ( |van den Berg et aT||2009| l. Our observations of the LW using the IMACS 
camera on the Magellan 6.5 m telescope show one or two candidate counterparts within the error circle of six 
sources, and all show relatively high X-ray-to-optical flux ratios. The observed periods (~ 1.3 to 3.4 hours) 
and the X-ray luminosities (io 31 - 8-32 - 9 erg s" 1 at 8 kpc) of the 10 periodic sources, combined with the lack 
of bright optical counterparts and thus high X-ray-to-optical flux ratios, suggest that they are likely accreting 
binaries, in particular, magnetic cataclysmic variables (MCVs). All of the 10 sources exhibit a relatively hard 
X-ray spectrum (r < 2 for a power law model) and X-ray spectra of at least five show an extinction larger 
than the field average expected from the interstellar medium in the region. The latter implies some intrinsic 
absorption in the system, which is also a typical sign of MCVs. The discovery of these periodic X-ray sources 
in the LW further supports the current view that MCVs constitute the majority of low luminosity hard X-ray 
sources (~ 10 30-33 erg s -1 ) found in the Bulge. The period distribution of these sources resembles those of 
polars, whereas the relatively hard spectra suggest that they could be intermediate polars (IPs). These p uzzling 
properties, w hich are also shared by some of the periodic X-ray sources found in the Sgr A* field ( |Muno 
et al.||20~03b[ l, can be explained by unusual polars with buried magnetic fields or a rare sub-class of MCVs, 
nearly synchronous MCVs. These unusual MCVs may provide important clues in the evolutionary path of 
MCVs from IPs to polars. While the Chandra X-ray band appears relatively more sensitive in discovering 
these unusual MCVs, the completeness simulation indicates >40% of the hard X-ray sources in the LW are 
periodic. Therefore, this discovery provides a first direct evidence of a large MCV population in the Bulge. 
Subject headings: Galaxy: bulge — X-ray: binaries — cataclysmic variables 



1. INTRODUCTION 



') on Galactic scales beyond the local 
Several ongoing campaigns, including 



The high sensitivity and superb spatial resolution of Chan 
dra enabled population study of low-luminosity X-ray sources 
(L x ~ 10 30 " 34 erg s 
solar neighborhood, 
our o wn Chandra M ulti-wavelength Plane (ChaMPlane) sur- 
vey ([Grindlay et al. 2005| aim to improve the census of the 
Galactic low-luminosity X-ray sources. The Galactic Bulge, 
in particular, has been of great interest. More than 3000 
discrete X-ray sources h ave been discovered in th e 17'xl7' 
region around Sgr A* (Mun cTet al. 1120 03a, 2009, hereafter 
M03a, M09). In ChaMPlane, we study the Bulge through our 
dedicated surveys of low-extinction bulge regions ("Windows 
survey") and a latitudinal strip around the Galactic center 
("Bulge Latitude Survey"), where thousands of X-ray sources 



2011 Hong et 



ha ve been also discovered ( |Grindlay et al 
aL"1|2009b| hereafter H09b). Multi-wavelength follow-up of 
these X-ray sources has been ongoing for source classifica- 
tion through optical/infrared imaging and spectrosc opy (e.g. 
IKoenig et al.|2008||van den Berg et al.|2006| 12009] hereafter 
B06,B09). 

H09b (see also M03a) have found that the X-ray sources in 
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the Galactic center region (GCR) show largely homogeneous 
X-ray properties (e.g. intrinsically hard X-ray spectra, pho- 
ton power law index T <1 for a power law model). The lack 
of bright IR (K<15) counterparts for the GCR X-ray sources 
indicates that HXMBs, once considered as a major constitute 
for Bulge X-ra y sources, cannot acc ount for more than 10% of 
the population ( |Laycock et al.|2003] hereafter L05). Currently 
the leading candidates that fit the observed properties are mag- 
netic cataclysmic variables (MCVs) (L05; M09; H09b). 

Recently we started a searc h for periodic ChaMPlane 
sources (e.g. Hong et aL]|2009a hereafter H09a), in part, to 
circumvent difficulties in source identification, which arise 
from large distances, high extinction (except for the Windows 
fields), and source confusion due to high stellar density in the 
Bulge. In this work, we report the discovery of 10 periodic 
X-ray sources from the 1 Ms Chandra exposure of the Lim- 
iting Window (LW), a low extinction region at 1 .4° south of 
the Galactic center (see also H09b). 

The low extinction Window fields, including the LW, 
provide a rare opportunity of studying the GCR source 
population without obscura tion from mol e cular clouds 
(H09a,H09b,B06,B09). |Revnivtsev et~all ( [20091 here- 
after R09) showed that the Galactic Ridge X-ray Emission 
(GRXE), the nature of which has been puzzling for decades, is 
mainly made up of discrete faint sources of known nature, pri- 
marily active binaries (ABs) or cataclysmic variables (CVs), 
based on the Chandra observations of the LW. However, the 
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TABLE 1 

Chandra ACIS-I OBSERVATIONS OF THE LW 



(JDS. ID 


Start Time 


K.A. 


Decl. 


KOli 


Oil 


Exposure 


epoch 




(UT: y/m/d h:m) 


(°) 


(°) 


(°) 


(ks) 


(ks) 


(duration) 


6362 


2005/08/19 16:15 


267.86875 


-29.58800 


273 


37.4 


38 




5934 


2005/08/22 08:16 


267.86875 


-29.58800 


273 


36.1 


41 












subtotal 


73.5 


79 


1 (3.13 d) 


6365 


2005/10/25 14:55 


267.86875 


-29.58800 


265 


20.4 


21 












subtotal 


20.4 


21 


2 (0.24 d) 


9505 


2008/05/07 15:29 


267.86375 


-29.58475 


82 


10.7 


11 




9855 


2008/05/08 05:00 


267.86375 


-29.58475 


82 


55.3 


57 












subtotal 


66.0 


68 


3 (1.22 d) 


9502 


2008/07/17 15:45 


267.86375 


-29.58475 


281 


161.7 


167 




9500 


2008/07/20 08:11 


267.86375 


-29.58475 


280 


158.6 


165 




9501 


2008/07/23 08:13 


267.86375 


-29.58475 


279 


130.0 


135 




9854 


2008/07/27 05:53 


267.86375 


-29.58475 


278 


22.4 


25 




9503 


2008/07/28 17:37 


267.86375 


-29.58475 


275 


99.5 


103 




9892 


2008/07/31 08:07 


267.86375 


-29.58475 


275 


64.8 


65 




9893 


2008/08/01 02:44 


267.86375 


-29.58475 


275 


41.4 


45 




9504 


2008/08/02 21:23 


267.86375 


-29.58475 


275 


122.8 


127 












subtotal 


801.2 


832 


4(17.7 d) 



"The selected Good Time Interval (GTI) that is based on the lack of high 
observations are grouped into four separate epochs (durations given in days), 



background fluctuations (< 3a). See Hong et al. 1 2005 i. The 
so that the total span of each epoch does not exceed a month. 



exact composition of the disc rete sources in the G RXE still 
remains unresolved (see also |Revnivtsev et al|201 We ex- 
plore the X-ray and optical properties of the 10 periodic X- 
ray sources in the LW and their implication for evolutionary 
models of MCVs and their connection to thousands of X-ray 
sources in the GCR. 

2. X-RAY OBSERVATION AND TIMING ANALYSIS 

2.1. Chandra Observation and Source Search 

The LW was observed for a total of 1 Ms exposure (100 ks 
in 2005 and 900 ks in 2008) with the Chandra ACIS-I instru- 
ment (H09b; R09). Table [T] lists the observational history of 
the field. The X-ray data were analyzed as a p art of our ongo- 
ing ChaMPlane survey (Grindlay et al. 2005 ) and the an alysis 
procedures are described in detail in |Hong et al.| ( [20Tj5| here- 
after H05; see also H09b). In summary, we stacked all of 
the 13 separate pointings and searched for discrete sources in 
the 0.3-2.5, 2.5-8. 0, and 0.3-8.0 keV ba nd images, using the 
wavdetect routine ( |Freeman et al.|[2002) >. The source lists of 
the three energy bands were then cross checked to produce a 
list of 1397 unique discrete sources based on the positional 
uncertainty of each source. Carefully designed aperture pho- 
tometry (H05) that takes into account source crowding was 
applied to each source in the 0.5-2.0, 2.0-8.0, and 0.5-8.0 
keV bands. The complete source list and their photometries 
will be presented elsewhere ( |Hong et al.||201l] ) along with 
comparison of the results from two popular source detection 
algorithms, wavdetect and wvdecom/Pj In this paper, we use 
the sources discovered by the wavdetect routine. Since we 
can identify periodic X-ray modulation only from relatively 
bright sources, the two search algorithms make no difference 
for the results of this paper. The catalog of the 319 sources 
and their photometry results from the initial 100 ks exposure 
can be found in H09b. 

2.2. X-ray Timing Analysis 
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The 1 Ms Chandra exposure of the LW can be divided into 
into four epochs so that the total span of each epoch does 
not exceed a month, which allows period search in a single 
ephemeris (Table fT). Of four epochs, the last epoch provides 
a long exposure (GTI: 801 ks) suitable for in-depth search of 
periodic modulation. For periodicity search, we selected the 
381 sources with background-subtracted net counts greater 
than 100 in the 0.3-8.0 keV band in Epoch 4 (Fig[]}. 

Photon arrival times of each source were Barycenter cor- 
rected to Barycentric Dynamical Time (TDB) by the axbary 
routine in the CI AO tool (ver 3.4). Then, we generated a 
background-subtracted light curve of each source i n 12.8 s 
bins, and we applied a Lomb-Scargle (LS) routine (Scargle 
1982) to the light curve to search for periodic X-ray modu- 
lation. We subtracted the background counts in the source 
aperture region using the events, with a proper scale, that fell 
in the background annulus region around the source, which 
excludes the source regions of neighbors (H05). The 12.8 s 
bin interval was chosen for efficient executions of the LS rou- 
tine over a long exposu re, a nd it only suppresses short periods 
below ~ 100 sec (see §6.1| and Fig. [7]), where there is already 
a concern due to the CCD readout time (~ 3.2 s). Random- 
ization of photon arrival times to compensate for fixed CCD 
readout cycle makes practically no difference for the search 
results due to the longer time bins used for the light curves 
and the long exposure consisting of multiple pointings, where 
multiple phases of CCD readout cycle produces an effect sim- 
ilar to randomizing arrival times to some degree. 

The search periods were selected successively, starting 
from the total duration (T) down to 1Q~ 4 T by a decrement 
of AP = P 2 /(2Ts f ), where we introduced an oversampling 
factor, Sf (s/=l means no oversampling), in order to sample 
the periodograms relatively smoothly over the entire period 
search range. We change Sf logarithmically from 1.0 at the 
shortest period (10~ 4 r) to 4.0 at the longest (T), and a rela- 
tive increase in the number of search periods due to the in- 
troduction of Sf is about 30%. In Epoch 4 where the total 
duration spans about 17.7 days with a combined total expo- 
sure of 832 ks, we searched 26499 periods from about 153 s 
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Total 


1397 


Counts>100 in Epoch 4 


381 


P FAP (Vor) < 1% 


17 


Pfap (Fix) < 1% 


14 


Periodic Sources 


10 



99% (Fix): 




99% (Vor) 



1000 

Net counts in Epoch 4 



10000 



FIG. 1 . — The distribution of the background subtracted counts in Epoch 4 
of the sources and the peak values of the periodograms by the LS method. Of 
17 sources with Pfap (Var) < 1% (net count dependent False Alarm Proba- 
bility; see j ]2.3| for its definition) or P con f (Var) > 99% (red solid), we identify 
10 periodic sources (see SRI. The periodicities of 5 sources (open circles) are 
false, originating from eitner the dithering motion of the instrument or non- 
periodic X-ray variabilities (e.g. flares). We also exclude two other sources 
('+' and 'x') because of the marginal nature of their detections (see 



to 17.7 days. In this analysis, we further limit the search range 
between 153 s and 10 hr, which contains 19911 independent 
periods out of 26247 trial periods. This range covers most of 
the spin and orbital periods of MCVs. We found that apparent 
periodicities of > 10 hrs indicated by the periodograms are 
false alarms, usually caused by short or long term X-ray vari- 
abilities (e.g. flares). For selected periodic sources (§2.3), we 
also extended the period search range down to 20 s in attempt 
to find any secondary modulation at periods shorter than 153 
s, although the bin size (12.8 s) of the light curves may sup- 
press detection of short periods below ~ 100 s (see Fig. IT!. 

For candidate periodic sour ces, we al so applied an Epoch 
Folding (EF) method ( Leahy et al.|1 983) to refine the periods 
and compare the significance of the primary periods with their 
harmonics. For the EF method, we generated a background- 
subtracted folded light curve in 15 bins for a given search pe- 
riod and calculated the x 2 value of the folded light curve with 
respect to the assumed constant count rate of no periodicity. 
For each source, we applied the EF method at 1000 periods 
within 2 a of the primary period found by the LS routine, and 
the same around the half, 2nd and 3rd harmonics of the pri- 
mary period. 



2.3. Periodic X-ray Source Selection 

In order to select valid periodic sources, we have esti- 
mated a False Alarm Probability (Pfap) or a confidence level 
(P con t = 1-Pfap)- For a given power density (X) in the LS 
routine, Pfap = 1 - (1 - e ) NsNp , where Np is the number of 
independent trial periods (19911) and N$ is the number of 
the sources. To validate Pfap estimated from the above trial 
statistics, we have also performed ~ 100,000 simulations by 
randomizing photon arrival times. These simulations allow 
reliable estimates of Pfap down to ~ 0.1% for Ns < 10 or 
~ 1% for Ns < 100, etc., where the simulation results match 
with those from the trial statistics. For a larger value of Ns or 
a smaller value of Pfap, more simulations are required. In the 
following, we calculate Pfap according to the trial statistics. 

The search is usually conducted for the brightest sources 
first, since they have a better chance of periodicity detection. 



Depending on the search limit of net count (N net ), which de- 
termines Ns, Pfap changes. Although the LS routine is known 
to be more sensitive to low count sources than a few other 
conventiona l Fourier analyses (e.g. see Chapter 6. 1 in |Bret-| 
|thorst|1988) , its performance eventually drops out as the net 
count decreases to zero. Extending the search limit by includ- 
ing low count sources with no chance of periodicity detection 
can arbitrarily increase Pfap, which can make valid periodic 
modulation appear insignificant. On the other hand, since 
search performance depends on various factors (e.g. back- 
ground counts) besides net counts, it is desirable to extend the 
search to a relatively low net count to be conservative. There- 
fore, we set the search limit at a relatively low value (N net = 
100), and we use two types of the Pfap estimates: Pfap (Var), 
where N s = N s (>N nel ), and P FA p (Fix), where N s = N s (N net 
>100) = 381. The former provides an initial cut for candidate 
periodic sources, and the latter provides a list of sources with 
a higher confidence. 

Modulation amplitude is also a good metric for significance 
of periodicity. In order to describe diverse pulse profiles, we 
use the rms amplitude (Arms) given by 
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where N and B are the total and background number of events 
in the source aperture region respectively, and tj is the ar- 
rival time, and <fi{t) is the phase at time t th at would be ex- 
pected for a given period ( Mu no et al.|200 3b hereafter M03b 
and the references therein^ For easy comparison with the 
literature, we also define another modulation amplitude as 
A M = 1 -Pmin/Pmax, where P min (P max ) is the minimum (max- 
imum) count rate of the folded light curve. Finally, for syn- 
thetic light curve simulations, we use sinusoidal variations, 
where the light curve is described as 1 +Aq sin </>(f), where 
< Aq < 1 and </>(f) is the phase of photon arrival time t. For 
sinusoidal modulations, Aq = v2Arms, = 2 Ao/(l +Arj), 
and Arms < Aq < Am- Since both P m ; n and P max depend on 
the bin size of folded light curves, we calculate A M from Arms 
using the above relations. 

We also use simulations to calculate a detection probability 
of periodicity (Pdet)- For each candidate periodic source, we 
generate 1000 synthetic light curves to estimate P^t within 
about 1% accuracy. Each synthetic light curve is consistent 
with a sinusoidal variation of the measured period and modu- 
lation amplitude of the source, while accounting for the GTIs 
and Barycentric shifts of the CCD readout time of the real 
data. Then we perform the LS search algorithm to see how 
often we detect the same periodicity with Pfap (Fix) < 1%. 
These simulations enable com pleteness correction of period- 
icity detection (see also §6.2| > and provide another validity 
measure of the detection in addition to Pfap- 

Finally, we also perfor med the Bay esian Block (BB) search 
for long term variability (Scargle 1998 ). In order to apply the 
BB search to the data set of a long duration (3 yr) with long 
exposure gaps (>1 yr), we eliminate exposure gaps longer 
than 20 ks, and generate a semi-contiguous series of photon 
arrival times. 

The simulations for Pfap and Pd e t along with period search 
require significant computational power. In order to handle 
the computational burden efficiently at low cost, we have uti- 
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TABLE 2 
Periodic sources in the LW 





1*91 






C51 


\S>) 


\ 1 ) 




(9) 


noi 


ID 


Snurpp "Mump 


ConritQ fFnnph Zt*l 


Period 


r FAP 


^RMS 


P i , 

* det 


Offset 


Signi. 


Notes 


LWP 


CXOPS 


(0 3-8 keVI 


(s) 






\ ,v ) 


\ ) 


ITfirmnnip*; 

1 1 HI 1 1 H 1 1 1 1 ^ n 




1 


J175151.2-293310 


1418(41) 


10342(5) 


19.6 


26(3) 


99.8 


5.6 






2 


J175123.5-293755 


323(21) 


5131(5) 


11.8 


59(5) 


99.9 


2.6 






3 


J175129.1-292924 


296(22) 


7448(10) 


6.8 


65(6) 


99.1 


6.1 






4 


J175131.6-292956 


2899(56) 


8536(14) 


6.1 


14(2) 


83.3 


5.6 




BB:4 


5 


J175133.9-292754 


419(29) 


6342(7) 


6.0 


43(4) 


67.5 


7.7 






6 


J175118.7-293811 


269(19) 


4729(2) 


4.0 


50(5) 


80.1 


3.3 






7 


J175122.7-293436 


f 119(13) 


12076(91) 


2.6 


60(8) 


25.7 


1.5 


2 


Edge: 352 ks 


8 


J 175 147. 4-2942 15 


307(28) 


4890(6) 


2.4 


40(4) 


9.6 


7.9 






9 


J175133.6-293313 


359(21) 


6597(10) 


2.1 


39(5) 


66.1 


2.6 


3 


PO:0.91 


10 


J1751 19.4-293659 


377(22) 


* 5262(2) 


1.4 


30(5) 


13.9 


2.4 


2 


2nd harmonic 



(1) An abbreviated X-ray source ID. (2) The Chandra source name. (3) The background subtracted net counts in the 0.3 -8 ke V band in Epoch 
4 (832 ks exposure). T The counts for LWP 7 is from the last 352 ks GTI, which is free of near CCD edge events (see i 5.7 1. See Table|3]for 
the total net counts of the 1 Ms exposure. 



(4) The modulation periods are refined by the EF rou tine a round the significant periods found by 
the LS routin e. *F or LWP 10, the second harmonic (5262 s) is considered the real period (see 5 5.10 1. (5) The False Alarm Probability for 
/Vi=381 (see j j2.3| l. (6) The rms modulation amplitude based on Eq. 1. (7) The detection probability of periodicity based on simulations using 
1000 syntheticTIght curves for each source. (8) The offset from the aimpoint of Obs. ID 5934 (R. A.: 17h 51m 28.50s, DEC: -29° 35'16.80"), 
(9) The most significant harmonic according to the E F method if the primary period is not. In LWP 10, the second harmonic is significantly 
more prominent than the primary primary (see j j5.10| l, whereas in t he other two sources, the differences are marginal. (10) BB: the number of 
independent Bayesian Blocks in the light curve if not one (see j |5.4} . Edge: the source falls near a CCD edge in some pointings and GTI free of 
such pointings is shown. PO: the 95% PSF overlaps with neighbors. The number indicates the fractional radius of the non-overlapping region 
(H05). 



lized a GPU-based desktop supercomputer equipped with a 
C1060 Tesla GPU (240 cores) from Nvida^] A simple code 
conversion using an IDL0GPU library (g/jM//^ has boosted 
the speed of some routines by a factor of 10 by comparison to 
the regular IDL routines^] 

3. X-RAY ANALYSIS RESULTS 

Fig. [T] illustrates the search results and the periodic source 
count. For initial screening, we select 17 sources that exhibit 
a periodicity with Pfap (Var) less than 1%. We exclude five 
sources with apparent periodicity due to either the dithering 
motion of the telescope (707, 1000 s and their immediate 
(sub-) harmonics), or non-periodic variabilities (e.g. flares). 
In addition, we exclude J175055.5-292948 (period: 5331 s) 
since its fdet is only 0.1% ('+' in Fig. [TJ. For sources with 
Pfap (Fix) > 1%, we only consider them periodic if they show 
another significant i ndicat or of periodicity (e.g. the results of 
the EF method, see ^10]l, and so J175103.9-293430 (period: 
3821 s) is also excluded fx' in the figure). The final list of 
the periodic sources in the LW contains 10 sources (9 with 
Pfap (Fix) < 1%). The spatial distribution of the 10 periodic 
sources does not show any obvious sign of clustering or pat- 
tern. 

Table [2]lists the basic results of period search and photom- 
etry of the 10 periodic sources in the LW. For easy reference, 
in addition to the X-ray source name (starting with CXOPS, 
see H05), we assign an abbreviated version of the name start- 
ing with prefix 'LWP' indicating 'Limiting Window Periodic 
sources'. The number IDs in the short name are assigned in 
the reverse order of the periodogram value at the observed pe- 
riods (i.e. the likelihood of true periodicity). Note B09 list the 
sources with prefix 'LW'. 



Figs. [2] and [3] show the periodograms, 
arves, tne X-ray spectra with model fits, ar 



the folded light 

curves, trie X-ray spectra with model fits, and the compressed 
light curves without long exposure gaps of the 10 periodic 
sources. The periodograms are based on the LS method, 
which was applied to the binned light curves of Epoch 4. The 
folded light curves are also generated from the data obtained 
in Epoch 4 and are very similar to the ones from the full 1 Ms 
exposure (not shown). The periodogram and the folded light 
curve of LWP 7 are from the last 352 ks GTI in 8. Id when the 
source was observed sufficiently distant of any CCD edges. 

The X-ray spectrum of each source is generated using all of 
the available data, and each spectral bin contains more than 20 
background-subtracted net counts. The red solid lines show 
the results using a simple power model fit, and for LWP 1, a 
spectral fit using a power law plus an iron line at 6.7 keV is 
also shown (blue). In the compressed light curves, the long 
exposure gaps (> 20 ks), are marked by the vertical lines, and 
the yellow lines indicate the boundaries of four epochs. The 
red horizontal lines show the average count rate and the green 
lines show the effective area, which varies from pointing to 
pointing. The blue solid lines for LWP 4 show the count rates 
of the BBs. 

Table [3] summarizes the X-ray spectral properties including 
estimates of the unabsorbed X-ray fluxes, based on spectral 
model fits. Note the net counts in Table |3] are from the full 
1 Ms exposure, whereas the net counts in Table [2] are from 
Epoch 4 except for LWP 7, which shows the net count from 
the last 352 ks GTI. We fit the spectra of all the sources using 
a simple power law model, and for LWP 1, we also tried a 



power law plus an iron line (see 5 5. 1 
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7 The complete code optimization is expected to boost the speed by another 
factor of 10 according to the examples given by gpulib. 



4. OPTICAL OBSERVATION AND PROPERTIES 

We observed the LW field with the HST and Magellan tele- 
scope. Five of the 10 sources were in the HST/ ACS fields 
(B09), and three of them were reported as possible candidates 
for accreting binaries or CVs. Table H] summarizes the com- 
bined optical and X-ray properties of the candidate optical 
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Periodogrom (LS) 
60 



Folded light curve (cts/ks) 



Energy spectrum (cts/s/keV & Sigma) 




" LWP 5 

P: 6341.54 s 
" (refined by EF) 


99% (Fix): 








99S (Vor): 





LWP 1 P: 10342-30 s 

A RMS : 26% (refined by EF) 



802ks (GTI) from 54664.6600 MJD (To) 



LWP 2 P: 5131.14 s 

A RUF; : 59% (refined by EF) 




802ks (GTI) from 54664.6600 MJD (To) 



LWP 3 P: 7448.40 s 

A RWS : 65% (refined by EF) 
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FIG. 2. — Periodograms, folded light curves, energy spectra and compressed light curves of the periodic X-ray sources (LWP 1-5) in the LW. The periodograms 
are based on the LS method, and the 99% confidence levels for N$ = 381 (Fix) and N$ = N$(> N n a) (Var) are shown in (red) horizontal lines (see j |2.3fr . The 
folded light curves are drawn with To = MJD 54664.6600, and the (red) horizontal lines show the average rates. In the energy spectra, the (red) lines show the 
results of simple power law model fits. In the case of LWP 1, a model fit using a power law plus an iron line is shown in blue. In the compressed light curves, the 
(red) horizontal lines represent the average rate, the (green) steps show the effective area of each pointing, and the vertical lines indicate the epoch boundaries 
(yellow) and exposure gaps (grey, > 20 ks). For LWP 4, the Bayesian Blocks are shown in blue lines. 



counterparts of the Chandra sources. For easy reference, the 
short names used by B09 are noted for the three sources in 
Table|4](e.g. LWP 2 = LW 25). 

4.1. HST/ACS Data and Analysis 

We observed with HST the inner area of the ACIS field 
with a 2 x2 mosaic of the Wide Field Camera (WFC) on the 
ACS on 2005 August 19. A single WFC pointing images a 
3.4' x 3.4' field with ~ 0.05" pixels using two CCD detectors 
separated by a 2.5" gap. Exposures were taken through the 
F435W (#435 ), F625W (R 625 , similar to Sloan r), and F658N 
(Ha) filters. Each tile of the mosaic was observed with the 
same exposure sequence 4 x492 s in F435W, 168 s + 2 x 167 
s in F625W, and 4 x496 s + 4 x492 s in F658N. No dither- 
ing was applied to fill in the WFC chip gap. Photometry is 
performed using a stellar-photometry package, DOLPHOT, a 
modified version of the HST phot package to do photometry 
on HST/WFPC2 images ( [Dolphin ||2000] >. See B09 for the 
details. 



4.2. Magellan/IMACS Data and Analysis 

On 2007 May 8, we observed the LW field and two other 
Window fields (Stanek's and Baade's) with the Inamori Mag- 
ellan Areal Camera and Spectrograph (IMACS) on the 6.5 m 
Magellan (Baade) telescope at Las Campanas, Chile. With 
seeing ~ 0.8- 1.2" FWHM, we obtained a dithered set of 5 
pointings in the f/4 configuration (15.4' field, 0.11'Vpixel) to 
cover an 18' x 18' region of the LW. This provided a total ex- 
posure time of 600, 300, 180, & 180 s in Bessell-B, V, R, & 
CTIO-7 filters over the Chandra field respectively. 

We processed the images using standard IRAF tasks, and 
calibrated the astrometry using the 2MASS catalog as a ref- 
erence. The astrometric residuals on each CCD frame were 
~ 0.2". We reprojected and stacked the images using the 
SWARFrjutility. All frames were normalized to ADU/second 
units anacombined using weight-maps constructed from flat- 
fields and bad pixel masks. The initial source search and 



http : / /www . astromatic . net / software /swarp 
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photometry were performed on the stacked images using 
DAOPHOT. See also H09a. 



4.3. Optical Matches 

Both the HST and IMACS source lis ts are bo r esigh ted to 
the Chandra sources as described in |Zhao et al.| ( |2005] ): the 
boresight correction is less than 0.1" in both R. A. and Dec- 
lination. In search of optical counterparts and their optical 
properties, for the five sources in the HST/ 'ACS fields, we use 
the HST observations. For the other five, we use the Mag- 
ellan/IMACS observations, where three sources show one 
or two candidate counterparts. Considering the high stellar 
density in the region, there is no guarantee of these Magel- 
lan/IMACS sources being the true counterparts (e.g. LWP 9 
has six candidate counterparts seen in the HST/ ACS field, but 
only one in the Magellan/IMACS image). 

The X-ray-to-optical flux ratios, Log(fjf/F^), in Tableplare 
calculated for the 0.3-8 keV band vs. the R magnitude. Both 
the observed and unabsorbed flux ratios are given. For unab- 
sorbed flux ratios, the R magnitud e is de-redde ned based on 
the Ap625\v map given by |Revnivtsev et al ( 2010 1. For sources 
with multiple candidate counterparts, a range of the R mag- 



nitude and the flux ratio is given, covering all the candidate 
counterparts. In the case of LWP 2 (= LW 25) and LWP 6 (= 
LW 8), shown is the R magnitude of the most likely candidate 
with an unusually blue color (B09). In the case of LWP 8, 
only one of the candidate counterparts has a measurable R 
magnitude. For sources with no detectable candidate coun- 
terparts, a lower limit for the R magnitude is given, based on 
the minimum value of optical neighbors within a 30" radius 
of the X-ray source position. 

The logarithmic flux ratios or their limits of these periodic 
sources are relatively high: > -1 for 7 sources and > -2.5 for 
the rest. The intrinsic X-ray-to-optical flux ratios of accreting 
binaries or AGN are usually significantly higher than those of 
coronal sources (see §4. 1 .2 in B09 and the references therein). 
CVs have the intrinsic flux ratio between -2.5 and +0.5 in the 
above energy bands, although a few outliers of active binaries 
or dMe stars have the flux ratios as high as -1 or -0.5. There- 
fore, the flux ratio results in Table [4] are consistent with those 
of CVs. 

Note that, except for a few cases with an outstanding blue 
counterpart, most of the flux ratio values in Tableware in fact 
likely lower limits, considering the high stellar density in the 
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TABLE 3 

X-RAY SPECTRAL PROPERTIES OF THE PERIODIC SOURCES IN THE LW BASED ON SIMPLE POWER MODEL FITS 



C1 1 


\^) 






(■>) 


W \i ) 




(8) 






OUUICC 














Unabsorbed Flux 


1 n 111 iTn~\ci1"\f 
J_,UllllIlUalLy 


ID 


Counts ( total) 


E50 


N H22 


1 




U.J— I 


2-8 


0.3-8 keV 


U.j— KeV 


LWP 


(0.3-8 keV) 


(keV) 






(keV) 




/1A-15 -2 

(10 erg cm 


s ) 


(10 x erg s" 1 ) 


1 


1775(45) 


3.86(6) 


1.4(2) 


0.4(1) 


0.86/74 


6.3(5) 


56(2) 


63(2) 


30.9-32.7 








1.7(1) 


0.63(2) 


0.5(1) 0.74/72 


8.3(7) 


55(1) 


64(2) 


30.9 - 32.7 


2 


393(23) 


2.9(1) 


2.2(5) 


1.9(3) 


1.07/14 


9(1) 


8.9(6) 


2KD 


30.4 - 32.2 


3 


361(24) 


3.7(2) 


1.4(9) 


0.8(4) 


0.74/13 


1.4(4) 


10.2(7) 


12.0(8) 


30.2 - 32.0 


4 


3217(59) 


3.36(4) 


1.5(1) 


1.12(8) 


0.94/123 


24(1) 


79(2) 


107(2) 


31.1-32.9 


5 


473(30) 


3.4(2) 


0.8(5) 


0.8(3) 


0.91/18 


2.5(4) 


12.7(9) 


16(1) 


30.3-32.1 
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3.2(1) 


1.3(6) 


1.2(4) 


0.35/11 


2.6(4) 


7.0(5) 


10.0(6) 
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0.32(6) 


7.8(6) 


7.6(5) 
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361(30) 


3.0(2) 
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10.3(9) 
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2.9(1) 
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0.8(3) 
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10 
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9.7(6) 


13.1(7) 
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(1) An abbreviated X-ray source ID. (2) The background-subtracted net counts in the 0.3-8 keV band from the full 1 Ms exposure. (3) The 
median energy of the photons in the 0.3-8 keV band. (4) & (5) An estimate of the power law index (r) and the absorption (^22) from spectral 
model fits, assuming a simple power law model for the X-ray spectrum in the 0.3-8 keV band. (6) The estimated equivalent width of the iron 
line from a spectral fit of a power law model with an iron line at 6.7 keV. (7) Reduced \ 2 and Degree of Freedom (DoF) of spectral model fits. 
(8) An estimate of the unabsorbed flux. (9) The luminosity range in the 0.3-8 keV band for source distance at 1-8 kpc. 



TABLE 4 

Candidate optical counterparts of the periodic sources in the LW 



(1) 


(2) 


(3) 




(4) 




(5) 


(6) 




(7) 


Source 
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B-R 


F x (0.3-8 keV) 


Log(F x /F R ) 


No. of Candidates 


AF625W 


LWP 


(mag) 




(observed) 


(unabsorbed) 


(observed) 


(unabsorbed) 


HST 


Magellan 
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1 


>22.3 




5KD 


64(2) 


>1.4 


>-0.1 






3.9(1.5') 


2 


22.3 .. 24.4 




8.0(5) 


21(1) 


0.6 .. 1.4 


-0.4 .. 0.4 


2 (LW25) 




3.6 (0.1') 


3 


>22.2 




9.2(6) 


12.0(8) 


>0.6 


>-1.0 






4.2 (1.9') 


4 


23.1(4) 




75(1) 


107(2) 


1.9(2) 


0.4(2) 




1 


4.1 (1.4') 


5 


20.5(2)* 




13.0(8) 


16(1) 


0.08(8)* 


-1.47(8)* 




1 


4.1 (3.4') 


6 


22.7(1) 


2.9 


6.8(4) 


10.0(6) 


0.68(5) 


-0.72(5) 


5 (LW8) 


1 


3.9 (0.3') 


7 


18.8 




7.6(5) 


7.6(5) 


-0.83(3) 


-2.47(3) 


1 


1 


4.1 (0.5') 


8 


21.4(2) 




8.1(7) 


10.3(9) 


0.23(9) 


-1.02(9) 




2 


3.4 (3.5') 


9 


22.4 .. 25.7 




8.6(5) 


9.3(5) 


0.7 .. 2.0 


-0.7 .. 0.6 


6 (LW19) 


1 


3.6 (0.2') 


10 


22.6 .. 24.6 




9.3(5) 


13.1(7) 


0.8 .. 1.6 


-0.6 .. 0.2 
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(1) An abbreviated X-ray source ID. (2) R magnitude, based on the HST observation if available: otherwise, based on the Magellan Images. 
For sources with multiple counterparts, the range of those counterparts is given. For sources with no valid counterparts, a lower limit is given, 
based on the minimum magnitude of neighboring optical sources within a 30" radius from the X-ray source position. '*' uses the / magnitude 
instead of R. (3) B (F435W) - R (F625W) from the HST observation. (4) The X-ray flux in 10~ 15 erg cnT 2 s~\ from column (7) in Table|] (5) 
Log(F x /F R ) = hog(Fx)+R/2.5 + 5J6. For unabsorbed flux ratios, the R magnitudes are de-reddened by Ap625w- (6) The number of candidate 
counterparts found in the Advanced Camera for Surveys (ACS) images of the HST and the Magellan MOSAIC Imag es. Every source in the 
ACS FoV has more than one counterpart and for sources reported by B09, their source IDs are given. (7) Af625w from[Revnivtse v~et al| ( |2010"| > 
and the offset of the source from where Af625w is sampled. 



region and the possibility of the true counterpart being fainter 
and undetected. LWP 5 and 8 have additional uncertainties 
in their estimates due to the variation in the interstellar ab- 
sorption across the region: for these sources, the Ap62sw val- 
ues were sampled about 3 or 4' away from the sources. For 
LWP 7 and 9, the absorption in the X-ray spectra is estimated 
less than what is expected in the field, based on A F 625w- For 
LWP 3, 5 and 8, two estimates are consistent, and for the rest, 
the former is larger than the latter. 

5. SOURCE PROPERTIES 

In this section, we discuss some of the unique properties or 
analysis caveats of each source. 

5.1. LWP 1: CXOPS J 175 151. 2-2933 10 

The periodogram reveals significant periodic modulation at 
the primary period and its half. The folded light curve shows 



a clear eclipse, but due to the lack of other features in the 
folded light curve other than the eclipse, the system can be 
either a polar or an IP. The long period suggests the modula- 
tion is likely due to the orbital motion (e.g. only two IPs with 
> 10 4 s spin period in Fig. [6^, see §6.1 for the details). If 
so, the mass of the companion is estimated to be ~ 0.42 M 
according to Eq. 2.89 in i Warner (1995). In addition, from 
the eclipse duration A<j) ~ 0.09, the mass ratio q > 0.29 and 
the inclination angle i is somewhere between 60° and 74° for 
w hite dwarf (W D) mass My ~ 0.1 - 1.4M (Eq. 2.92 & 2.93 
Wamer|1995) . 



m 



The X-ray spectrum exhibits a feature consistent with an 
iron line at 6.7 keV (and perhaps 6.4 keV as well). The spec- 
tral fit with a power law plus a 6.7 keV line reduces the %?, 
by 14%, compared to a simple power law model (Table |3f. 
In ord er to estimate the significance of the line (see Protassov 
2002} , we generated 1000 synthetic spectra, each of which is 
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consistent with a simple power law model (T = 0.4 and Nmi 
= 1 .4). We then fit each spectrum with a power law plus the 
iron line, and count how many cases reduce the xt by more 
than 14%. The results were none, indicating the significance 
of the line is 99.9% or higher. The hard X-ray spectrum with 
an iron emission line suggests the source is likely an IP (see 



more about the source type in §6.1 



5.2. LWP2: CXOPS J 17 5123. 5-2937 55 

This source was reported as a possible accreting binary 
(LW 25) by B09. Two optical sources in the HST/ACS images 
are found within the error circle of the X-ray source position, 
and neither of them stands out with any unusual colors. The 
folded light curve shows a possibility of a long eclipse start- 
ing at phase ~ 0.5. LWP 2 and 6 are located in the region with 
an excess of seemingly diffuse soft X-ray background, which 
perhaps indicates the lowest extinction region of the field. In 
the case of LWP 2, it shows a significantly larger absorption 
in the X-ray spectrum than the field average, indicating an in- 
trinsic absorption in the system. 

5.3. LWP 3: CXOPS J 17 5 129. 1-292924 

This source exhibits a mild anti-correlation between the 
count rate and the absorption in the X-ray spectr um accord- 
ing to a phase-resolved quantile diagram (Fig. |4} ( jHong et aT] 
|2004| H09a). This anti-correlation is consistent with a picture 
that the observed X-ray modulation is caused by the variation 
of intrinsic absorption in the system, similar to the IP discov- 
ered in BW (H09a). Alternatively if the X-ray modulation 
originates from an eclipse or an obscuration of the hot spot 
or the emission region due to the spin or orbital motion, then 
the absorption variation may not be expected to be strongly 
correlated with the rate change. The folded light curve of the 
source shows a hint of an eclipse (phase ~ 0.45), which, if 
true, implies the system is sychronized (i.e. polar), consider- 
ing the synchronized primary modulation in the folded light 



5.4. LWP 4: CXOPS J 1751 31. 6-292956 

This is the brightest source of the 10 periodic sources. The 
X-ray spectrum does not show any sign of iron lines. The 
BB search indicated four independent blocks in the long term 
light curve, but since the source fell near a CCD edge in Epoch 
3, only two blocks can be credited to be independent. In sum- 
mary, the observed X-ray flux in 2008 is higher by a factor of 
- 1.8 than that in 2005. 

5.5. LWP 5: CXOPS J 17 5 133. 9-2927 54 

This source and LWP 8 were observed at relatively large 
off-axis angles. 

5.6. LWP 6: CXOPS J1751 18.7-29381 1 

This source was reported as a candidate CV (LW 8) by B09, 
based on the high X-ray-to-optical flux ratio and Ha excess 
(see Fig. 4 in B09). The folded light curve in Fig. [2] shows 
a hint of an eclipse at phase ~ 0.2. Our routine analysis of 
phase-resolved quantiles at the primary period and the second 
harmonic revealed an intriguing result as illustrated in Fig. [5] 
At each period, we show two versions of folded light curves 
using sliding phase windows: one with fixed-width phase bins 
(10%, grey), and the other with variable-width phase bins but 
fixed net counts in each phase bin (45 counts, black). The 



former is better suited to identify a sudden drop in the count 
rate such as an eclipse, and the latter is better suited to reveal 
a sudden increase. In addition, the latter is better suited for 
phase-resolved spectral quantile analysis since every phase 
bin contains enough events to allow a reliable estimate of en- 
ergy quantiles. The folded light curves in fixed-width phase 
bins show an eclipse-like feature both at the primary period 
(phase ~ 0.2 and 0.6) and at the second harmonic (~ 0.25). 
The ingress and egress of the eclipse are sharper at the second 
harmonic than the primary period. A narrow eclipse, synchro- 
nized with the primary modulation, suggests the system is a 
polar. 

What is interesting in this source is a spectral change cor- 
related with the phase. At the primary period, the first quar- 
ter of the phase shows an intrinsically harder X-ray spectrum 
(T < 1) than the rest (T > 1), which implies two different 
emission regions or mechanisms are present (Fig.[5£). On the 
other hand, at the second harmonic, there appears to be a dra- 
matic change in the absorption between the first (Nmi ^ 2) 
and second halves (Nmi ;$ 1) of the phase (Fig.[5}l & f). The 
change in the absorption between two peaks of the folded light 
curve would suggest that two magnetic poles of the system are 
visible in turn, and the X-ray emission from one of the two 
undergoes through more material (likely an accretion curtain, 
trailing to a pole and extended from about half of the accre 
tion disk or ring) before reaching us. Fig. 5 in Evans & Hel 



lier ( 2007 ) illustrates a possible viewing geometry for such 



a system: unlike their exmaples, the soft blackbody compo- 
nent (<C 1 keV) of LWP 6 is likely always invisible due to the 
interstellar absorption, but the geometry allowing for a phase- 
dependent variation of intrinsic absorption due to an accretion 
stream or curtain can apply to LWP 6. This picture, if true, 
suggests that the second harmonic (9457 s) is a real period and 
the observed primary period (4729 s) is a sub-harmonic, even 
though the LS method did not find the second harmonic sig- 
nificant. Further observation is required to determine which 
period represents the true orbital and spin geometry. In addi- 
tion, a new diagnosis may be required to quantitatively eval- 
uate the signifi cance of various types of spectral changes (see 
|Connors|20TT] l. 

5.7. LWP 7: CXOPS J 17 5122.7 -293436 

During three pointings out of 8 total in Epoch 4, this source 
fell near a CCD edge (see the effective area in Fig. [3]), which 
would discredit the observed periodicity, but the relatively 
clean data set (the last 352 ks GTI) free of near-CCD-edge 
events also exhibits a significant periodicity at 12076 s. The 
net count of LWP 7 in the last 352 ks GTI is only 119, 
which is similar to that of J 175055.5-292948 ('+' in Fig.[T} in 
Epoch 4, but Pd et of LWP 7 is estimated substantially higher 
than J175055. 5-292948 because of the relatively lower back- 
ground counts in the source aperture region of LWP 7. There- 
fore, we consider the observed periodicity of LWP 7 valid. 

5.8. LWP 8: CXOPS J17 5147 .4-294215 

The folded light curve of the source shows an eclipse-like 
feature at phase ~ 0.6. 

5.9. LWP 9: CXOPS J175 133.6-2933 13 

LWP 9 was reported as a potential accreting binary (LW 19) 
by B09. The aperture source region (95% PSF) of these 
sources mildly overlaps with that of a neighbor. A clean data 
set relatively free of contamination from the neighbor (see 
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FIG. 4. — Phase-resolved quantile analysis of LWP 3 using the last 801 ks exposure: (a) the folded light curves, (b) the rate vs. the quartile ratios, and (c) the 
phase-resolved quantile diagram. The folded light curves are generated using sliding windows of fixed-width phase bins (grey) and variable width phase bins 
with fixed net counts (black: example data points with bin sizes are shown in phase 0—1) and the error bars of the latter are color-coded by phase for easy 
comparison. The quartile ratios is mostly proportional to Log(A^22) f° r LWP 3. 
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FIG. 5. — Phase-resolved quantile analysis for LWP 6: the left panels for the primary period (4731 s) and the right panels for the second harmonics (9461 s). 
(a,b) the folded light curves, (c,d) the phase-resolved estimates of using the power law models in the quantile diagram and (e,f) the phase-resolved quantile 
diagram. For the 9461 s period (d), the first half of the phase (0.0-0.5) shows a higher absorption (N^22 ~ 2) than the second half (0.5-1 .0) (N^22 ~ 1), 
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H05 for aperture choice) exhibits the same periodic modu- 
lation, so that its periodicity is considered valid. 

5.10. LWP10: CXOPSJ1751 19.4-293659 

According to the EF analysis for LWP 7, 9 and 10, a higher 
harmonic exhibits more significant periodicity than the pri- 
mary period found by the LS method. While the differences 
in the significance between the primary periods and their har- 
monics are marginal for two sources, in the case of LWP 10 
the second harmonic appears substantially more significant. 
The x 2 of the 15-bin folded light curve with respect to the 
constant rate is 133 at the second harmonic vs. 89.5 at the pri- 
mary periocj^] Therefore we consider the second harmonic to 
be the real modulation period for LWP 10. The high \- value 
at the second harmonic is partially due to an eclipse-like fea- 
ture in the folded light curve at phase ~ 0.7, which suggests 
the system is a polar. 

6. DISCUSSION 

In this section, we investigate the most probable source 
types for the periodic sources found in the LW and their impli- 
cations, based on the statistical distribution of their properties. 
We also explore the hidden population of periodic sources in 
the Bulge X-ray sources in the GCR through completeness 
simulations for periodicity detection. 

6.1. Unusual MCVs? 

The observed source properties such as the X-ray luminos- 
ity range (~ 10 30-33 erg s" 1 for distance of 1 - 8 kpc), the 
relatively hard X-ray spectra (T < 2), the period range (1.3 - 
3.4 hr), and the relatively high X-ray-to-optical flux ratios, all 
indicate that these sources are typical MCVs. But the collec- 
tion of these properties does not appear to fit well with most 
common types of MCVs as explained below. 

MCVs can be largely divided into two groups, IPs and po- 
lars, depending on the relative strength of the magnetic fields. 
Traditionally, polars are synchronized or nearly synchronous 
(PJP ~ 0.98 - 1.02), whereas IPs are not (PJP < 0.1). Figffl 
shows the spin vs. orbital period distribution of the MCVs 
in the Ritter & Kolb (RK) catalog (ver 7.15) ( |Ritter & KoTb 
2003) >. In Fig. [6^, the observed periods of the periodic X-ray 
sources found in the LW and Sgr A* fields are shown as or- 
bital periods along PJP a = 2 or 4 for easy comparison. 

In the case of the periodic sources in the LW, the observed 
period distribution (red in Fig. [SJ)) resembles those of polars 
(green) better than either the spin (black) or orbital (blue) peri- 
ods of IPs. Some IPs do have spin or orbital periods at around 
an hour to three hours, but the majority of the spin (or orbital) 
periods are shorter (or longer), whereas the majority of the pe- 
riods of polars are in the same range as those of the periodic 
sources in the LW. 

In order to find out if the above result is due to a period- 
dependent selection bias in the periodicity search routines, 
we have conducted a set of simulations using synthetic light 
curves with various net counts (100 to 1000), modulation am- 
plitudes (A = 10% to 100%) and periods (~ 20 to 1.1x10 s 
s). For each combination of parameters, we generate 500 syn- 
thetic light curves, which allows < 2% accuracy in measure- 
ment of the detection probability. As in the simulations in 

9 For comparison, the expected \ 2 of the 15-bin folded light curve is 70.9 
at Pfap = l%forJYs=381 and Np= 199 11, although the EF search was con- 
ducted only for the 10 sources at 4000 periods around the likely true modu- 
lation periods. 
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each synthetic light curve is generated to 
properly reflect the GTIs and Barycentric shifts of the CCD 
readout time of the real data. 

The simulation results show there is no significant selection 
bias in the range of 150 sec to 10 hr (Fig. 17). This implies 
that the resemblance to the polar period distribution is not due 
to any selection effect in the search algor ithmsj^] or if these 
sources are IPs, they indeed belong to a statistically differ- 
ent population from the typical IPs in the RK catalog. In ad- 
dition to the observed eclipses (or eclipse-like features) syn- 
chronized with the primary modulation of the pulsed profiles 
(e.g. LWP 3, 6, 8 or 10), missing secondary periods from all 
the sources in Table [2] indirectly supports the systems being 
polars - synchronized systems, although non-detection does 
not guarantee the absence of the secondary periods. 

Fig [8] illustrates a wide range of modulation amplitudes of 
the 10 periodic sources in the LW. For comparison, we also 
show some of literature-selected IPs and polars. If the X- 
ray emission originates from a small spot (e.g. polar cap) on 
the WD surface, the modulation due to the spinning motion 
of the compact object is expected to exhibit a larger ampli- 
tude change than that from the orbital motion. For instance, 
a sample of IPs in Fig. [8] selected from the literature, show 
a slightly higher average value of the modulation amplitude 
at the spin periods (black closed circles) than at the orbital 
periods (black open circles). However, the modulation distri- 
butions of these IPs and polars in Fig. [8] are likely selection 
biased: e.g. 17 polars in Fig. [8] are mostly eclipsing systems, 
which are likely preferentially found in periodicity searches 
due to the large modulation amplitude. Note for polars in 
Fig. [8] Am is simply based on R m i„ and R max of the published 
light curves, where by definition, Am=\ for eclipsing systems, 
whereas for the periodic sources in the LW, we calculate Am 
from A rms (e.g. for LWP 1, Am=0.53 instead of 1.0). In sum- 
mary, it is not easy to link the distribution of the modulation 
amplitudes to a particular type of MCVs due to the lack of sys- 
tematic survey. However, even in Fig.|8]with selected samples 
of IPs and polars, the period distribution of the 10 periodic 
sources strongly favors polars. 

Unlike the period distribution, the relatively hard X-ray 
spectra of the periodic sources in the LW imply that these sys- 
tems are likely IPs. The quantile diagrams in Fig. [9] illustrate 
the overall spectral hardness of the periodic sources (red) in 
comparison with the rest of the X-ray sources (black dots, net 
counts > 50) in the LW. Fig. [9h and b overlay a set of sim- 
ple power law and APEC model grids respectively over the 
same data points. The eight periodic sources found in the Sgr 
A* field are also shown in (green) crosses (M03b). All of the 
periodic sources show an intrinsically hard X-ray spectrum, 
similarly to the periodic sources in the Sgr A* field. 

IPs tend to show harder X-ray spectra, which are associ- 
ated with higher acc retion rates and weaker magnetic fields 
(Ramsay & Cropper 2004). For instance, in the case of po- 
lars, the X-ray spectra are well described by a blackbody 
component with kT < 60 eV and a two-temperature thermal 
plasma compon ent with kT\ = 0.7-0.9 keV and kT2 = 3-5 
keV (e.g. Rams ay et al.|20 04b), whereas the X-ray spectra of 
IPs show a blackbody component with kT > 60 eV and a one 
or two-te mperature th ermal plasma component with kT > 10 
keV (e.g. |Anzolin|2008[ ). Due to the interstellar absorption in 
the LW field (A^h22 ~ 0.7), the blackbody component is usu- 

10 There is, however, a selectio n bia s towards sources with high modula- 
tion amplitudes, as expected. See 36.21 
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FIG. 6. — (a) The spin (P s ) vs. orbital period (P ) distribution of MCVs from the RK catalog (ver 7.15). The vertical yellow lines mark the period gap, and the 
diagonal lines represent Pa/P = 1, 0.6, 0.3 and 0.1. The purple diamonds show eight asynchronous polars (APs) (three are on top of each other), and there are six 
near synchronous IPs with P s /P > 0.3. The observed periods of the periodic sources in the LW (red) and the Sgr A* field (blue) are marked as orbital periods, 
along P s = 2 P or 4 P D for clarity (those less than an hour are likely the spin periods), (b) The period distribution of the LW sources (red) is compared to the 
distribution of the spin (black) and orbital (blue) periods of IPs and the periods (green) of polars in the RK catalog. The periodic sources (red) in the LW fill in 
the period gap, indicating these are likely MCVs. The periods of the LW sources are distributed more closely to those of polars (green) than the spin (black) or 
orbital (blue) periods of IPs. 
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FIG. 7. — Examples of detection (Pfap (Fix) < 1%) probabilities of peri- 
odicity as a function of period by simulations. The plot for "All Combined" 
includes the following cases: for net = 1000, An = 0.2 .. 0.5, for net = 500, 
A = 0.3 .. 0.7, for net = 200, A = 0.5 .. 1.0, for net = 150, An = 0.6 .. 1.0, 
and for net=100, An = 1.0. The detection probability show no bias within the 
period range of the interest between 150 s to 10 hr (two vertical grey lines). 



ally undetectable, but the spectral distinction of the plasma 
components between polars and IPs remains detectable in the 
Chandra X-ray band: in a quantile diagram, polars would lie 
in the upper-left section of kT ~ 4- 10 keV line, and IPs in the 
lower-right section (Fig. [9]). Although the above description 
of the X-ray spectra of IPs and polars is over-simplified and 
witho ut a systematic survey (e.g. LWP 6 is likely a polar, see 
^5.6[ ), it is generally accepted that IPs exhibit a harder spec- 
trum. For instance, a recent survey of hard X-ray sources (> 
15 keV) conducted by Swift IBM and INTEGRAL/IBIS shows 
that the composi tion of MCVs i n the hard X-ray band are pre- 
dominantly IPs (Scarin gT|2010| l. 

The period range and the relative hard X-ray spectra of the 
LW periodic sources also resemble those of some of the pe- 
riodic sources found in the Sgr A* field (green data points in 
Fig. [9] and blue points in Fig. |6^) (M03b). Of eight periodic 
sources in the Sgr A* field, four are in the same period range 
as the LW periodic sources with a large modulation amplitude 
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FIG. 8. — The X-ray modulation (Am) vs. the period distribution of 
literature- selected IPs (black closed circles for spin periods and black open 
circles for orbital periods), polars (black 'x's), and the periodic sources in 
the LW (red circles), BW (purple diamond) and t he Sgr A* field (green di- 
amonds). [References] Orbital periods of IPs: Parker, Norton, & Mukai 
(2005) . Spin periods of IPs: GK Per, HT Cam EX Hya, RX J 1548.2, AO Psc, 
V 1223 Sgr, AE Aq r -[Evan s & Hellief1 {2005) , V4 05 Aur -[Ev ans & Hellier 
1 2004 1, FO Aqu -|Evans et al.H2OTT 7 PTTTJem -|james et al.H2002) , IGk 
7131794 -[Butters et al.H2009l DU Col -|de Martino et al.N1006|rWX Pyx - 
Schlege f|j2005f , XV ARl -fcali nas & Schlegel ( 2004 li TTe^ -[Vr ielmann 
|et al.|(20(j4}7TV Col - |Rana et alH2004), V1062 Tau -|Hellier, Be ardmore 
& Mukai ( 2002a}, 1WG A J1958~2 -|Norton et al.H2002j , YY bra - Szkody 
et al.N 2002 1, V709 Cas - Norto nl^!!!^}. Po IarsTV1309 Ori -|Schwarz 
|et al.||2Tj 05 ), EK UMa - Beuermann et al. ( 2009 1, HU Aqr - |Schwarz eTaT 
TOT9TV2301 Oph - Ramsay & Cropper|(2o07} , SPSS JO 15543 -[S chmidt 
POTB) , E P Dra -|Ramsay et al.H2004c| , OYTarjWheatley cfWestPTOT 
UP L eo -|Ramsay et al.[12Uulf , RX J 1 8 46.9 -|Schwarz et al.H2TJiJ21 , V407 
Vul - Marsh & Steeghs ( 2002 1, CE Gru - |kamsay &"Tro pper 1 2002 1 V1432 
- |Rana et al|ft 005 1, V34 7"Pav, GG Leo, E U UMa- Ramsay et al. ( 2004a I, 
RX J 1002-19 - |Ramsay & Cropper|j20oi) ) 



(■Arms >40% or Am >72%), and all of them show relatively 
hard X-ray spectra. One of them was identified as a fore- 
ground polar, based on the light curve, and it lies at the line 
of kT = 10 keV and the rest lie in the kT > 10 keV section. 
In the case of the Bulge X-ray sources in the Sgr A* field, 
the heavy interstellar absorption (~ 6x 10 22 cm" 2 ) may com- 



12 



6 keV 




0.5 0.6 0.7 

log 10 (E 50 /0.3 keV)/log 10 (8/0.3) 



0.5 0.6 0.7 0.8 

log, (E 50 /0.3 keV)/log, (8/0.3) 



FIG. 9. — Quantile diagrams (H04, H()9a) of the periodic X-ray sources in the LW (red) and the Sgr A* field (green). The same data points are plotted over a set 
of power model grids (a) or APEC model grids (b). The (black) dots in (a) are the discrete X-ray sources in the LW with net counts greater than 50 in the 0.3-8 
keV band. The model grids are based on the response function at the aimpoint, which can be partially responsible for a small discrepancy between the spectral fit 
and quantile analysis in the parameter estimates for some sources (e.g. in LWP 9, Nh22 = 0-2(3) from the spectral fit vs. 0.7(2) from the quantile analysis). 



plicate identification of the intrinsic X-ray spectra, but for the 
periodic sources in the LW (~ 7xl0 21 cm" 2 ), there is little 
doubt that most of them exhibit intrinsically hard X-ray spec- 
tra. 

An obvious question, then, is, what is the nature of the peri- 
odic X-ray sources found in the LW and the GCR, whose pe- 
riod distribution resembles polars' but X-ray spectra resemble 
IPs'? 

First, one can consider these are a rare type of polars with 
unusually hard X-ray spectra. Although uncertain, it is spec- 
ulated that the origin of harder X-ray spectra in IPs relative to 
polars is related to the weaker magnetic field strength and the 
subsequently deeper penetrati on of the accretion stream into 
the WD surface. For instance, Cummming (2002) suggested 
the relatively high accretion rate in IPs effectively buries the 
WD magnetic field, making them appear less magnetic, which 
in turns helps maint aining the hig h accretion rate (e.g. V407 
Vul: see Marsh & Steeghs 2002 ). This idea nicely ties the 
dichotomy of the magnetic field strength between IPs and po- 
lars. Under this picture, one can imagine some unusual evo- 
lutionary scenarios from IPs to polars, where polars survive 
with a high accretion rate that can bury the magnetic fields, 
allowing a harder X-ray spectrum. If most of these periodic 
sources are located near the Galactic center (8 kpc) as ex- 
pected from the high absorption in their X-ray spectra and the 
high stellar density of the Bulge, the X-ray luminosities of 
these sources are estimated at the high end of the MCV range 
(> 10 32 erg s -1 , see Table[3]), which is also consistent with the 
above picture. 

Second, a rare subclass of MCVs, nearly synchronous 
MCVs (ns-MCVs), perhaps also meet both of the observed 
properties - the period distribution and the relative hard X- 
ray spectra of the periodic sources in the LW. One can divide 
ns-MCVs in two subgroups - nearly-synchronous IPs (ns-IPs) 
and asynchronous polars (APs). Both subgroups may ex- 
hibit similar X-ray properties, but probably represent different 
stages in the evolutionary path of MCVs. 

First, APs, consisting originally of just four systems, which 
recently extended to eight according to the latest RK catalog 
(ver 7.15), are traditionally considered as polars that are tem- 
porarily out of synchronization due to a recent Nova activity, 
which has altered their magnetic locking, giving PJP Q ^0.98- 



1 .02. Interestingly it is speculated that APs exhibit a harder 
spectrum than normal polars, but with the similar periods, as 
marked with purple diamonds in Fig. [6^. For instance, two of 
seven APs as opposed to two of 92 normal polars are found in 
the hard X-ra y survey using I NTEGRAL/IBIS, Swift/EKT and 
Suzaku/HXD ( |Scaringi|20T0l ). 

Second, there are increasingly more IPs found near at syn- 
chronization (PJP a > 0.3). Starting with EX Hya, the list 
increases to six according to the RK catalog. Their orbital pe- 
riods are predominantly around 1.5 hr except for V697 Sco 
with a 4 .5 hr orbital period . According to the evolutionary 
model of |Norton et al.| ( |2008| hereafter N08), IPs start out with 
PJPo < 0.1 and as the systems evolve through magnetic lock, 
the orbital periods decrease and the spin periods increase, i.e. 
PJPo approaches 1. Therefore, the orbital periods of ns-IPs 
will be clustered around or below the period gap near the end 
of the evolution, resembling the period distribution of polars 
more closely than that of usual, unsynchronized IPs. 

As with polars exhibiting unusally hard X-ray spectra, the 
presence of ns-MCVs is very intriguing in terms of the evolu- 
tionary models of MCVs, challenging the conventional view 
of IPs with PJP a < 0.1 an d polars ~0.98 - 1.0 2. For instance, 
Paloma or RX J0524+42 ( |Pineault et al.|1987| ), recently iden- 
tified as an AP, shows PJP ~ 0.93, and its relatively large de- 
synchronization (~ 7%) compared to the conventional APs (< 
2%) suggests that this system might represent the missing link 
of the evolutionary path between IPs and polars ( Schwarz et 



al. 2007[ ), rather than polars being out of synchronization tem- 
porarily due to recent Novae activities. As shown in Fig. [6^, 
the gap (~ 0.3 - 0.95) in the period ratio between IPs and po- 
lars are now more or less bridged by the recent discoveries of 
many ns-MCVs. It is speculated that some of these ns-MCVs 



transit between IPs and polars (e.g. V1025 Cen, see Hellier. 



Wynn & Buckley]: 

According to N08, if the period ratio exceeds 0.6, t he only 
stable equili brium is at synchronization (PJP =1). |Schwarz| 
|et al.| ( [2007] ) suggest that the probability of finding ns-MCVs 
is very low, considering relatively short timescale for syn- 
chronization (< 1 Myr) compared to the lifetime of a CV (~ 
100 Myr) (see also N08). Therefore, if many of the periodic 
sources in the LW are indeed ns-MCVs, it imposes another 
constraint on the evolutionary model or suggests an unusual 
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environment of the Galactic Bulge capable of harboring many 
such rare systems. The similar statement can be made for po- 
lars with unusally hard X-ray spectra. 

The relative composition of source types in MCVs are 
highly biased, depending on search wavelengths. In the RK 
catalog, where most of the discoveries are based on opti- 
cal/UV or longer wavelength bands, the relative ratio of IPs 
vs. polars are close to 1 (e.g. 83 IPs vs. 109 p olars in F ig.|6l>, 
whereas a hard X-ray survey (> 15 keV) in (Scaringi 20IT5) 
revealed 37 IPs and only 2 polars (both are APs). Therefore, 
if many of the periodic sources in the LW are indeed unusual 
polars or ns-MCVs, it indicates the Chandra X-ray band is 
well tuned for discovery of these rare MCVs. 

Finally we note that some of the observed X-ray properties 
of the periodic X-ray sources in the LW are shared by a group 
of quiescent Low Mass X-ray Binaries (qLMXBs). For in- 
stance, Swift J1353. 5-0127, recently observed in an outburst 
as a black hole (BH) trans ient, may have a rela tively short or- 
bital period (~ 2 hr, see |Casares et al. 2011|l and exhibit a 



hard X-ray spectrum (Krimm, Kennea & Holland 201 1 ). The 
observed over-abundance of X-ray transi ents within 1 pc of 
Sgr A* (cf. 4 out of 7 within 20 pc, see |Muno et al.H2005) 
implies that a large number of qLMXBs ma y be present as 
dark stellar remnants within 1 pc of Sgr A* ( |Schodel et al. 
2007). However, CVs (and MCVs) are expected dominantly 
more abundant than these BH transients, and the periodic X- 
ray sources in the LW do not exhibit any strong outbursts in 
the 1 Ms exposure spanned over 3 years (cf. out of 7 X-ray 
transients in the Sgr A* field, 3 or 4 sources have been ob- 
served in outbursts in each year, see Degenaar & Wijnands 



2009 2010) it is reasonable to think that the large fraction of 



the periodic sources in the LW are in fact MCVs. 

6.2. Periodic source content 
in the GCR X-ray source population 

In this section, we estimate the total number of periodic 
X-ray sources in the LW through completeness simulations 
for periodicity detection. We have generated 500 synthetic 
light curves for a given set of X-ray modulation parameters in- 
cluding modulation period and amplitude, net count and back- 
ground count. Then we run the same detection algorithm to 
see how often the synthetic light curves are detected as peri- 
odic. This calculates P^ et for the given set of parameters. Here 
"detection" is simply defined as Pfap (Fix) < 1 %FH The sim- 
ulation is designed to match the data set observed in Epoch 
4 with the same GTI gaps, and the simulated time tags are 
generated according to the Barycenter corrected CCD read- 
out time cycles. The initial set of simulations indicates that 
the detection results do no t dep end on the given period range 
(100 s to 10 hr) (Fig.[7]in §|6T}, so we fix the period at 5432.1 
sec and varied the rest of the parameters. The simulated mod- 
ulation amplitude (Ao) ranges from 10 to 100% in 10% incre- 
ments. The simulated net counts of sources are 50, 100, 150, 
200, 500, and 1000, 3000, and the background counts are 30, 
50, 100, 200, 500, and 1000. 

Using these simulation results, which cover most of the pa- 
rameter space for the sources found in the LW, we estimate 
the periodicity detection probability (Pdet) of the 381 sources 
with net counts > 100 in Epoch 4. Since we do not know 
a priori the distribution of modulation amplitude and the ob- 
served periodic sources show a wide range of the modulation 

1 1 The result in this section is consistent with that acquired with "detection" 
defined as P FAP (Var) < 1%. 



amplitude (Ao), we randomly assign a modulation amplitude 
(Ao) to each source, assuming a uniform distribution from 0. 1 
to 1. We repeat the interpolations 10000 times (i.e. 10000 set 
of simulations), each with random assignments of modulation 
amplitude to cover the full modulation range for every source. 

Fig. 10 shows an example result, using 203 sources with net 
counts > 100 and £50 > 2.5 keV (hard X-ray sources, see also 
case 3 in Table Bl) under the assumption that these sources are 
all periodic with"! a uniform distribution of the modulation am- 
plitude. Fig.flOfr shows the expected number of sources with 
detectable periodicity in the 2-D phase space of net counts 
vs. modulation depth for hard X-ray sources. The absolute 
number in each cell of the phase space depends on the cell 
size, so what matters here is a relative variation from cell to 
cell and the total number. The result shows that we expect to 
detect periodicity of 21(3) sources if all of the 203 sources are 
periodic with a uniform distribution of modulation amplitude. 
Fig.~ 

The source distribution 



10 3 & c show the same result (differential and accu 
mulated) as a function of net count 



with detectable periodicity (red) indicates that the detection 
probability for sources with net counts less than 150 is very 
lowp] Since we only detected 9 periodic sources with Pfap 
(Fix)< 1%, this result indicates there should be about 87^2 
periodic sources, which is 43^% of the 203 sources. How- 
ever, since the detection probability drops significantly at net 
counts below 200, we cannot make a reliable assessment of 
the faint periodic source population. If we limit our detection 
probability at 1%, the dotted line in Fig.[T0p shows the effec- 
tive population of the sources we can explore. Accordingly 
the actual estimates of the periodic population can be as low 
as the dotted green line in Fig. 10 ;. Therefore, in the worst 
case scenario where most of the sources with net counts below 
200 are not periodic, the estimated periodic source population 
is 50^7°, which is about 24^% of the hard X-ray sources. 

Fig. [Toll shows the accumulated percentage of the periodic 
sources in the LW as a function of net count. Fig.[l0}l shows 
a gradual increase of the fractional periodic sources as net 
counts decrease to 200 from 500. This observed variation of 
the fractional periodic sources with net counts is consistent 
with the fact that the X -ray fluxes of the m ajority of MCVs are 



below 10 erg s ( Heinke et al.|2008| ) since the observed 
periodic sources are at the high end of the X-ray luminosity 
distribution if they are in the Bulge. It also implies the true 
percentage of the periodic sources with net counts less than 
200 is likely higher than the value (43%) at net count ~ 200. 

Table [5] summarizes the simulation results that estimate 
the total number of periodic sources including unidentified 
ones among the sources found in the LW. We repeated the 
analysis for four sub-sets of the sources, which cover sam- 
ples only with non-CCD-edge sources and hard X-ray sources 
(£50 > 2.5 keV). In this table, we show the results for sources 
with net counts greater than 200, where we have a sufficiently 
high detection probability of periodicity given the number of 
sources. The case using non-CCD-edge sources enables an 
estimate of the systematic error in our selection procedure of 
periodic sources since the observed periodicity of the sources 
that fall near a CCD edge can be falsely discredited. The case 
using only hard X-ray sources should be a better representa- 
tive of the Bulge sources since they exclude the majority of 
the foreground soft sources. 

12 This is an over-simplification since Pfe t also depends on the background 
counts. See LWP 7 in Tableland ^377] 
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FIG. 10. — Completeness study for periodicity detection for hard X-ray sources (£50 > 2.5 keV; see also case 3 in Table|5j. (a) The (relative) expected distri- 
bution of detectable periodic sources based on the completeness simulation under the assumption of all the sources being periodic with the uniform distribution 
of modulation amplitudes. The open circles show the 10 periodic sources in the LW. (b) The number of the sources with detectable periodicity (red) as a function 
of net count if all the sources (black) are periodic in comparison with the observed periodic sources (blue). The dotted line indicates the testable sample of the 
population if we limit the detection sensitivity at 1%. (c) The accumulated number of the detectable sources (red) if all the sources are periodic (black). Given 
the observed periodic sources (blue), we also show the expected total population of periodic sources (green). Below 200 counts, the detection sensitivity is too 
low to make a reliable assessment of the hidden population of periodic sources (dotted green lines). The shade for green and blue indicates the statistical errors 
of the simulation results and estimates, (d) The fractional percentage (solid line) of the periodic sources, which is the ratio of the green to black lines in (c). The 
dotted line shows a conservative limit below 200 counts due to lack of detection sensitivity, but the rapid rise of the percentage of peri odic sources as the count 
approaches from 500 to 200 indicates the trend will continue at low counts, where the majority of the Bulge MCV populations lie (see j |6.2) . 



TABLE 5 

Estimation of the total number of periodic sources in the LW by completeness simulations 



(1) (2) 


(3) 


(4) 


(5) (6) 


(7) 


(8) (9) 


Cases Source Selection 


Total 




Periodic X-ray Sources 




Alternative Estimates 


(net > 200 in Epoch 4) 


Sources 


Observed 


Detectable Estimated Total 


Percentage 


Estimated Total Percentage 


1 all 


153(12) 


9 


28 (-4, +4) 48 (-6, +8) 


32 (-5, +6) 


48(15) 32(10) 


2 non-edge 


113(11) 


8 


21 (-4, +4) 43 (-6, +9) 


38 (-7, +9) 


43(14) 38(12) 


3 £50 > 2.5 keV 


96(10) 


9 


21 (-3, +3) 42 (-6, +8) 


43 (-8, +10) 


42(12) 43(13) 


4 £50 > 2.5 keV, non-edge 


69(8) 


8 


16 (-3, +3) 35 (-6, +8) 


51 (-11, +14) 


35(11) 51(16) 



(1) Case number (2) Source selection criteria. We limit the sample to sources with > 100 net counts in Epoch 4. (3) The total number of 
sources that meet the source selection criteria. (4) The number of the observed periodic sources from Table [2] (5) The estimated number of 
the detectable periodic sources by simulations under the assumption that all the sources are periodic with a uniform distribution of modulation 
amplitude. (6) The estimated total number of the periodic sources in source selection: (3) x (4) / (5). (7) The percentage of periodic sources 
in the selection: (6)/(3). (8) An alternative (error) estimate of the total number of periodic sources and (9) their percentage using negative 
binomial distribution, where, for the observed periodic sources, N, and the non-detection probability, p, given by the simulations (in case 3, 
N = 9, l—p = 21/100), the estimated unidentified periodic sources and its variance are Np/(l-p) and Np/(l-p) 2 respectively. 



A simple ratio argument between the observed periodic 
sources and simulation results (column 6 & 7 in Table[5]) indi- 
cates that about 32-38% of the sources in the LW should be, 
in fact, periodic and the percentage increases to 43-5 1 % for 
the hard X-ray sources. The errors of these estimates are cal- 
culated from the variances of the total sample sizes (column 3) 
and the ranges of the detectable periodic sources (column 5), 



provided that at least the observed number of periodic sources 
(column 4) are present in each case. Alternatively one can es- 
timate the errors for the total number of the periodic sources 
according to negative binomial distribution (column 8 & 9 in 
TableBl), which provides the variance estimates of the uniden- 
tified (and thus total) number of periodic sources. 
The result indicates a large fraction of the Bulge X-ray 
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sources (> 30-40% of hard X-ray sources even with errors) 
found in the GCR should be periodic. This provides a first di- 
rect evidence for the presence of a large population of MCVs 
in the Bulge X-ray sources, which supports that MCVs con- 
stitute a large majority of the Bulge X-ray sources. 

7. SUMMARY 

We have found 10 periodic X-ray sources in the LW from 
the 1 MS exposure of the Chandra ACIS-I instrument. The 
overall period distribution, X-ray luminosity and spectral 
properties of these periodic sources and eight periodic sources 
found in the Sgr A* field (M03b) fit the general description 
of MCVs, supporting the argument that MCVs are the ma- 
jor constituent of the Bulge X-ray sources. However, when 
inspecting the details of these X-ray properties - period dis- 
tribution that resembles the polars' and the hard X-ray spectra 
that resemble the IPs', we discover that these periodic sources 
may fit into rare MCVs - polars with unusally hard X-ray 
spectra or nearly synchronous MCVs. The completeness tests 
indicate at least about 40% of the hard X-ray sources in the 
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LW are periodic, which is a first direct evidence of a large 
MCV population in the Bulge X-ray sources. 

Given the difficulty of identifying the nature of the Bulge 
X-ray sources in the GCR due to high stellar density, large 
distance and heavy obscuration, we recommend continuous 
X-ray monitoring of the GCR for future discoveries of period- 
icity and unusual variability. Considering the importance of 
these unusual MCVs in understanding the evolutionary path 
of MCVs and their potential connection to thousands of the 
Bulge X-ray sources, we also encourage systematic in-depth 
X-ray studies of the similar kind (e.g. there are only 13 ns- 
MCVs with PJP > 3 in the RK catalog) to establish a demo- 
graphic profile of their X-ray properties if any. 
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